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Structure  of  the  Compressible  Turbulent  Shear  Layer 


Dimitri  Papamoschou  * 

California  Institute  of  Technology,  Pasadena,  California 


Abstract 


Introduction 


The  large-scale  structure  of  the  turbulent  compressible 
shear  layer  is  investigated  in  a  two-stream  supersonic 
wind  tunnel  through  a  series  of  experiments.  Double¬ 
exposure  schlieren  photography  reveals  that  the  two 
convective  Mach  numbers,  corresponding  to  each  side 
of  the  shear  layer,  are  very  different,  one  sonic  or  super¬ 
sonic  and  the  other  low  subsonic.  This  contradicts  the 
current  isentropic  model  of  the  structure  which  pre¬ 
dicts  them  to  be  equal  or  very  close.  It  is  shown  that 
addition  of  shock-wave  effects  to  that  model  allows  for 
the  asymmetric  trends  observed  in  the  experiments. 
An  inclined  view  of  the  flow  provides  sketchy  informa¬ 
tion  about  the  spanwise  orientation  of  the  large-scale 
structure  and  does  not  reveal  any  pronounced  obliq¬ 
uity.  Attempts  to  enhance  mixing  by  modifying  the 
trailing  edge  were  unsuccessful. 


Nonmenclature 

a  speed  of  sound 

M  free  stream  Mach  number 

Mc  convective  Mach  numer  (unspecified  side) 

U  free  stream  velocity 

Uc  convective  velocity 

P  obliquity  angle 

7  specific-heat  ratio 

S'  growth  rate 

0  boundary-layer  momentum  thickness 

p  density 

Subscripts: 

1  high-speed  stream 

2  low-speed  stream 
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It  is  a  widely-accepted  premise  that  shear-layer  en¬ 
trainment  and  mixing  are  governed  by  the  instability 
of  the  turbulent  large-scale  structure.  The  visualiza¬ 
tion  of  such  structures  in  the  compressible  case  led 
Papamoschou  and  Roshko1'2  to  characterize  the  com¬ 
pressibility  of  the  flow  in  a  frame  of  reference  in  which 
the  structure  is  stationary.  Given  the  flow  conditions 
depicted  in  Fig.  1,  this  was  done  by  defining  the  con¬ 
vective  Mach  numbers  as  follows: 


(1) 


where  Uc  is  the  convective  velocity  of  the  structure, 
assumed  to  be  constant. 


= 


Ul-Ue 


a  i 


Me,  = 


Uc-Ui 

a2 


The  relation  between  Me ,  and  Mc,  was  obtained  by 
requiring  that  the  total  pressures  of  the  two  streams 
in  the  convective  frame  be  equal.  This  stems  from 
a  well-known  argument,  used  primarily  with  subsonic 
flows,  that  there  exists  a  stagnation  point  between  two 
structures  that  must  be  stable,  thus  the  pressures  at 
that  point  must  balance  (Fig.  2).  It  has  been  pro¬ 
posed  explicitly  or  implicitly  by  a  variety  of  investiga¬ 
tors  (Fliigel3,  Brown4,  Coles5,  Dimotakis6). 


The  above  argument  was  first  extended  to  compress¬ 
ible  flow  by  Bogdanoff7.  This  extension,  used  subse¬ 
quently  by  Papamoschou  and  Roshko,  assumes  that 
the  flow  comes  to  rest  at  the  stagnation  point  isenirop- 
ically,  hence  with  no  shock-wave  losses,  even  for  su¬ 
personic  convective  Mach  numbers.  For  streams  with 
equal  specific- heat  ratios  (71  =72)  one  simply  gets 


Mtl 


=  Me,  = 


A  V 

ai  +  02 


(2) 


where  A U  =  U\  -  U2.  For  the  cases  where  71  ^  72, 
Me,  and  Mc,  are  slightly  unequal  but  (2)  is  still  a  good 
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approximation.  Using  (1)  and  (2),  an  expression  for 
Uc  is  obtained: 


Uc  -Ui  1 

(/.-I/.-  1  +  VC 


(3) 


which  is  identical  to  that  obtained  by  the  subsonic 
argument.  The  experiments  of  Brown  and  Roshko8 
demonstrated  the  validity  of  (3)  for  incompressible, 
variable-density  flow. 

Since  MCl  and  A/c,  are  theoretically  equal  or 
very  close,  Papamoschou  and  Roshko  chose  MCl  as 
a  compressibility-effect  paramater  against  which  the 
shear-layer  growth  rate  was  correlated.  Their  result 
is  shown  on  Fig.  3,  where  the  growth  rate,  S' ,  is  nor¬ 
malized  by  its  incompressible  value,  6'0.  The  fact  that 
all  the  experimental  points,  corresponding  to  a  large 
range  of  free  stream  conditions,  collapsed  roughly  onto 
one  curve  demonstrated  the  applicability  of  of  the  con¬ 
vective  Mach  number  as  a  measure  of  shear-flow  com¬ 
pressibility. 

Nevertheless,  questions  remained  as  to  the  accuracy 
of  (2)  and  (3)  when  shock  waves  form  on  the  struc¬ 
ture,  since  their  effects  are  not  addressed  by  the  above 
model.  There  is  little  doubt  that  shock  waves  will  oc¬ 
cur  as  the  convective  Mach  number  exceeds  some  crit¬ 
ical  value.  This  is  the  basic  motivation  for  the  work 
described  here,  in  which  the  convective  Mach  num¬ 
bers  were  measured  experimentally.  Additionally,  ex¬ 
periments  related  to  the  spanwise  organization  of  the 
large-scale  structure  and  to  mixing  enhancement  were 
conducted  and  are  described  in  later  sections  of  this 
paper. 


Experiment 


illustrated  on  Fig.  4.  A  detailed  description  of  the  sys¬ 
tem  will  be  provided  in  a  different  publication10.  It 
retains  the  basic  elements  of  a  schlieren  system,  only 
now  there  are  two  slightly  off-axis  beams  produced  by 
two  different  spark  gaps.  By  means  of  the  prism  shown 
on  the  schematic,  the  source  focal  points  of  the  two 
beams  coincide  at  the  single  knife  edge  and  two  dis¬ 
tinct  images  appear  on  the  film.  The  exposure  time 
of  each  image  is  20  ns,  short  enough  to  capture  the 
details  of  the  flow. 

The  time  interval  between  spark  firings  (and  hence 
betwen  images),  At,  is  controllable.  During  that  time 
interval,  a  large-scale-structure  feature  moves  a  dis¬ 
tance  Ax.  The  convective  velocity  is  simply  Uc  = 
Ax /At.  The  convective  Mach  numbers  are  then  ob¬ 
tained  from  (1). 

The  measurement  error  for  Uc  is  roughly  ±5%, 
based  on  the  uncertainty  in  locating  a  given  feature 
from  one  exposure  to  the  next.  For  Mc ,  this  translates 
to  an  error  of  the  order  of  ±5 %/Mc.  Clearly,  for  Mc 
close  to  zero,  this  is  a  substantial  error.  But  for  Mc  of 
the  order  of  one,  the  error  is  similar  to  that  for  Uc. 

Fig.  5  shows  actual  measurements  of  Uc  from  two- 
spark  schlieren  photos  for  three  of  the  ten  Mach 
number-gas  combinations  explored  in  this  study.  Ar¬ 
rows  indicate  the  structure  features  on  which  the  mea¬ 
surements  are  based.  These  photos  were  taken  with 
the  knife  edge  parallel  to  the  flow  direction  and  span  a 
distance  from  50mm  to  140mm  from  the  trailing  edge. 
Other  knife  edge  orientations,  such  as  45-deg.  and  per¬ 
pendicular,  were  also  tried  and  did  not  produce  any 
difference  in  the  Uc  measurements  (see  Fig.  14).  For  all 
cases,  Ue  was  found  to  be  independent  of  streamwise 
position.  The  standard  deviation  of  the  measurements 
is  within  the  experimental  error  margin. 


Experiments  were  conducted  in  the  GALCIT  25mm 
x  57mm  x  200mm  supersonic  shear  layer  facility,  ex¬ 
tensively  described  in  previous  reports1,2,9.  Typical 
flow  static  pressures  were  0.1  atm.,  with  the  resulting 
unit  Reynolds  numbers  of  the  order  of  104  per  mm.  In 
the  region  of  the  measurements,  x/6\  exceeded  1000, 
where  is  the  trailing-edge  momentum  tickness  for 
the  high-speed  stream.  The  flow  is  thus  believed  to  be 
fully-developed  turbulent. 

Convective  velocity  measurements  were  taken  by 
means  of  a  two-spark  variant  of  the  schlieren  method, 
developed  by  the  author.  Its  operating  principle  is 


Results 

The  table  that  follows  summarizes  the  experimen¬ 
tal  conditions  and  convective  Mach  number  measure¬ 
ments.  Ten  cases  of  different  Mach  number-gas  com¬ 
binations  are  listed.  For  brevity  and  ease  of  refer¬ 
ence,  each  case  is  assigned  a  code  name  consisiting 
of  the  gases  and  Mach  numbers.  Letters  represent  the 
gases:  A  for  argon  (7=1.67),  H  for  helium  (7=1.67), 
N  for  nitrogen  (7=1.4),  and  S  for  sulfur  hexafluoride 
(SFe,  7=1.09).  Numbers  represent  the  Mach  numbers 
times  10.  For  example,  A32S03  is  the  case  with  argon 
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at  Mi=3.2  and  SFe  at  M?=O.Z.  The  faster  stream  is 
placed  first  in  the  name. 

The  cases  are  listed  in  order  of  increasing  .  As 
noted  previously,  this  parameter  is  very  close  or  equal 
to  the  theoretical  MCl  and  Mc„  so  it  replaces  them 
in  the  table.  The  last  column  contains  the  number  of 
data  points  ( N )  reflected  in  each  measurement. 


CASE 

'  Us,  ' 
Vi 

£2. 

Pi 

m . 

Ol  +  02 

MCl 

Mc, 

N 

A32N16 

0.94 

0.24 

0.07 

0.07 

0.07 

5 

N31N17 

0.75 

0.54 

0.33 

0.29 

0.36 

6 

N28A26 

0.75 

1.8 

0.38 

0.48 

0.26 

8 

H17N28 

0.50 

9.2 

0.65 

0.83 

0.10 

10 

A32A02 

0.13 

0.23 

0.90 

0.39 

1.14 

8 

S27S03 

0.13 

0.67 

1.06 

0.42 

1.61 

4 

H26N28 

0.42 

5.5 

1.09 

1.47 

0.10 

16 

H31N16 

0.30 

2.5 

1.38 

2.00 

0.32 

4 

A32S03 

0.08 

0.83 

1.55 

0.13 

3.15 

2 

N30S03 

0.06 

1.87 

1.71 

0.44 

3.67 

7 

The  experiments  revealed  that,  at  high  compressibil¬ 
ity  (case  H17N28  and  below),  Uc  closely  approaches  U\ 
or  U 2,  depending  on  the  test  case.  This  is  a  large  devi¬ 
ation  from  (3)  and  produces  the  surprising  convective 
Mach  number  trends  seen  on  the  table.  One  gets  a 
better  view  of  that  discrepancy  when  the  experimen¬ 
tal  Mc ,  and  MC}  are  plotted  versus  each  other.  This 
is  done  in  Fig.  6,  where  the  theoretical  values  are  also 
included.  Dashed  lines  connect  the  theoretical  and  ex¬ 
perimental  data  for  each  case.  The  slight  deviation  of 
the  theoretical  points  from  a  perfect  diagonal,  implied 
by  (2),  is  due  to  the  small  effect  of  different  7’s.  The 
difference  between  theory  and  experiment  is  striking: 
instead  of  being  equal  or  close,  Mc,  and  Mc ,  are  very 
different,  one  low  subsonic  and  the  other  sonic  or  su¬ 
personic. 

Equally  surprising  is  the  fact  that  the  side  of  the 
shear  layer  with  the  higher  Mc  is  not  the  same  from 
case  to  case.  In  some  cases  MCl  is  higher  while  in 
others  Mc,  is  higher.  There  is  no  consistency  of  these 
trends  based  on  quantities  that  are  frame-of-reference 
independent,  namely  on  **  and  on  For  example, 
cases  H31N16  and  N30S03  have  similar  ^  and  2a,  yet 
their  Me  distributions  are  different. 

However,  there  is  an  apparent  consistency  based 
on  the  free  stream  Mach  numbers:  in  supersonic- 
supersonic  combinations,  MCl  is  always  highest  (Ue  —* 
C/2)  while  in  supersonic-subsonic  combinations,  Me,  is 


always  highest  ( Uc  — • ►  U\ ).  This  baffling  trend  could  be 
a  coincidence  or  the  result  of  a  subtle  flow-field  change 
that  occurs  when  the  low-speed  stream  goes  from  sub¬ 
sonic  to  supersonic.  A  possible  reason  for  such  change 
is  given  in  the  discussion  section. 

Revised  Structure  Model 

It  is  evident  that  the  relationship  between  MCj  and 
Mc,,  derived  from  the  existing  large-scale  structure 
model  (Fig.  2),  is  grossly  inaccurate  for  high  compress¬ 
ibility  levels.  The  basis  of  that  model,  however,  is  still 
believed  to  be  sound:  in  the  convective  frame,  fluid 
engulfed  into  the  mixing  zone  has  to  come  to  rest  at 
some  point  between  two  structures,  and  the  pressures 
there  must  balance  for  that  point  to  be  stable.  The 
assumption  of  a  stationary  stagnation  point  could  be 
debatable.  Since  the  flow  is  turbulent,  the  stagnation 
point  could  be  shifting  around,  thus  adding  a  time  de¬ 
pendence  to  the  equilibrium  pressure  condition.  Judg¬ 
ing  from  the  two-spark  schlieren  photos,  however,  this 
effect  seems  small.  For  example,  in  the  photo  of  case 
H17N28  (Fig.  5)  the  structure  features  are  remarkably 
frozen  from  one  exposure  to  the  next,  although  they 
are  experiencing  a  shear  of  600  m/s. 

It  is  the  process  by  which  fluid  comes  to  rest  at 
the  stagnation  point  that  needs  critical  examination. 
The  existing  model  assumes  this  process  to  be  isen- 
tropic.  To  make  the  model  more  realistic  for  the  com¬ 
pressible  case,  the  dissipative  effects  of  shock  waves 
must  be  added.  A  shock  wave  formed  on  the  structure 
could  cause  considerable  total-pressure  drop  along  the 
streamline  leading  to  the  stagnation  point.  If  shock 
waves  appear  symmetrically  on  the  structure,  that  loss 
would  be  similar  on  both  sides  and  (2)  would  still  be 
a  good  approximation. 

There  is,  however,  potential  for  asymmetric  situ¬ 
ations.  A  “cat’s  eyes”  idealization  of  a  large-scale 
structure  sytem,  depicted  in  Fig.  7,  helps  visualize 
the  following  argument.  We  place  ourselves  in  the 
convective  frame  of  reference  and  assume  that  the 
two  free  streams  have  equal  static  pressures.  In  that 
frame,  suppose  that  the  Mach  number  of  the  upper 
free  stream,  Mc,,  is  supersonic.  As  in  flow  around  a 
cylinder  or  a  thick  airfoil,  the  Mach  number  near  the 
structure  is  higher  than  the  free  stream  Mach  number. 
Flow  along  the  top  edge  of  the  structure  turns  and  ac¬ 
celerates  through  a  Prandtl- Meyer  expansion  fan.  As 
a  result,  the  Mach  number  ahead  of  the  stagnation 
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point  is  substantially  higher  than  MCl.  Flow  comes 
to  rest  at  the  stagnation  point  by  means  of  a  nor¬ 
mal  shock  located  a  short  distance  upstream  of  that 
point.  A  normal  shock  at  that  high  Mach  number  is 
strong  enough  that  the  total  pressure  downstream  of 
the  shock,  i.e.,  the  pressure  experienced  by  the  stag¬ 
nation  point  on  that  side,  is  reduced  to  the  level  of 
the  static  pressure  in  the  free  streams.  As  an  exam¬ 
ple,  consider  a  case  with  MCl  =  1.0  and  7i=1.4.  For 
simplicity,  assume  that  the  Mach  number  at  the  top 
of  the  structure  equals  Afc, .  Suppose  that  the  shock 
occurs  at  a  flow  inclination  of  40  deg.  with  respect  to 
the  free  stream.  From  the  Prandtl-Meyer  function,  the 
Mach  number  before  the  shock  is  found  to  be  2.4,  suf¬ 
ficiently  high  that  the  total  pressure  after  the  shock 
equals  the  free  stream  static  pressure. 

In  the  convective  frame  again,  the  free  stream  Mach 
number  of  the  lower  side,  MC3,  must  be  such  that  the 
stagnation  point  is  pressure  balanced.  Two  situa¬ 
tions  are  now  possible-  a  symmetric  one,  already  noted 
above,  in  which  the  lower  flow  is  also  supersonic  with 
MC3  sb  MCl,  with  a  similarly-strong  shock  forming  on 
the  lower  side  of  the  structure;  an  asymmetric  one,  in 
which  the  lower  flow  is  subsonic  (Afc,  «  0),  with  the 
resulting  total  pressure  on  that  side  close  to  the  free 
stream  static  pressure.  Since  no  shocks  are  possible 
at  that  low  Mach  number,  the  total  pressure  remains 
virtually  constant  along  the  streamline  leading  to  the 
stagnation  point  and  balances  the  total  pressure  pro¬ 
duced  by  the  upper  side’s  shock.  Hence,  the  stagnation 
point  is  still  stable,  although  AfCl  and  MCi  are  very  dif¬ 
ferent.  These  two  situations  are  illustrated  in  Fig.  7, 
where  their  respective  total-pressure  distibutions  along 
the  stagnation  streamline  are  sketched. 

Obviously,  the  same  argument  can  be  considered 
with  MCo  given  supersonic,  in  which  case  one  gets 
AfCl  sb  MCi  (symmetric)  or  AfCl  ss  0  (asymmetric).  So 
a  total  of  three  situations  are  possible,  one  symmetric 
and  two  asymmetric. 

It  is  thus  seen  that  inclusion  of  shocks  in  the  system 
can  dramatically  alter  the  condition  for  pressure  bal¬ 
ance  of  the  stagnation  point.  Although  not  proposed 
as  the  only  possible  explanation,  the  revised  model  ra¬ 
tionalizes  the  basic  experimental  trends  produced  here. 
In  its  present  form,  however,  it  cannot  predict  which 
side  of  the  layer  will  have  the  higher  Me.  Also,  it  does 
not  explain  why  given  the  symmetric  and  asymmet¬ 
ric  possibilities  at  high  compressibility,  the  flow  here 
always  prefers  the  asymmetric  one. 


Significance  of 


It  is  important  to  develop  a  means  of  estimating  the 
likelihood  of  asymmetric  situations.  According  to  the 
revised  model,  flow  along  the  structure  boundary,  as 
seen  by  the  structure,  must  be  supersonic  for  such  sit¬ 
uation  to  be  possible.  Consequently,  the  free  stream 
convective  Mach  number,  Afe,  must  exceed  some  crit¬ 
ical  value,  M*.  As  "noted  earlier,  flow  in  the  vicinity 
of  the  structure  is  faster  than  in  the  free  stream,  so 
Af’  could  be  less  than  1.  A  reasonable  range  for  Af* 
is  0.5  <  Af *  <  1.  For  an  asymmetric  situation  to  be 
possible,  we  must  then  have  MCl  >  Af*  or  MCj  >  A/* . 


The  pertinent  question  to  ask  is  under  what  condi¬ 
tions  will  the  flow  have  no  choice  but  to  have  at  least 
one  Oi  the  Mc’s  greater  than  A/*.  To  help  answer  that, 
the  diagram  of  Fig  8  is  drawn.  It  shows  the  depen¬ 
dence  of  MCl  and  MC3  on  Uc  from  their  definitions  (2). 
Ue  is  bounded  by  the  free  stream  velocities,  U 2  and  U\. 
The  question  is  now  rephrased  as  follows:  regardless 
of  what  value  Uc  takes  in  the  interval  what 

criterion  ensures  that  MCl  or  MCi  or  both  will  exceed 
M*?  The  answer  is  obvious  from  the  diagram.  The 
Afc,  and  MCi  lines  form  a  boundary  whose  minimum 
occurs  at  MCl  =  MCl  =  ^+£7-  Thus  the  criterion 


A  U 

ax  -f  a2 


m; 


(4) 


guarantees  that,  no  matter  what  the  value  of  Uc  is,  at 
least  one  of  the  Mc’s  will  exceed  Af*  and  that  asym¬ 
metric  solutions  will  therefore  be  likely. 

It  is  evident  from  Fig.  8  that  asymmetric  solutions 
are  still  possible  when  a^fa  <  Af*.  In  that  regime, 
however,  the  flow  has  the  choice  of  a  subsonic  struc¬ 
ture,  that  is,  one  with  both  Mc’s  less  than  Af*  and  to 
which  the  isentropic  model  applies.  Given  the  sym¬ 
metric  choice  with  subsonic  flow  and  the  asymmetric 
choice  with  a  strong  shock  wave,  it  is  unlikely  that  the 
flow  will  prefer  the  highly-dissipative  asymmetic  one. 

Judging  from  the  results  of  this  study,  a  good  guess 
for  Af*  is  0.6.  The  experimentally-encountered  flow 
states  can  then  be  summarized  as  follows: 


<0.6 


Symmetric  situations 


>  0.6 

Oj  +aa 


Asymmetric  situations 
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The  parameter  originally  a  byproduct  of  the 

isentropic  model,  now  takes  a  significance  of  its  own. 
It  defines  the  possibilities  for  asymmetric  or  symmetric 
states,  and  does  so  in  a  way  largely  independent  of  a 
specific  large-scale  structure  model.  It  can  be  consid¬ 
ered  as  a  broad  measure  of  shear-flow  compressibility, 
since  its  value  does  not  depend  on  shear-layer  side. 

The  actual  physics  of  the  mixing,  however,  depend 
on  the  individual  convective  Mach  numbers.  In  partic¬ 
ular,  they  depend  strongly  on  the  higher  the  two  Mc’s, 
denoted  A/Cmar .  The  lower  Mc  is  always  low  subsonic, 
so  its  influence  is  weaker.  There  is  a  interesting  cor¬ 
relation  between  MCnmi  and  ,  shown  in  Fig.  9, 
where  a  monotonic  trend  is  observed. 

Structure  Obliquity 

Implicit  in  the  preceding  analysis  and  discussion  is 
the  assumption  that  the  spanwise  turbulent  structure 
is  two-dimensional.  It  is  possible,  however,  that  the 
structure  is  oblique,  with  a  direction  of  propagation 
inclined  at  an  angle  f3  to  the  mean  flow  direction  The 
compressibility  of  such  structure  is  governed  by  an  ef¬ 
fective  convective  Mach  number 

Mc.j,  =  MCC0S(3  (5) 

rather  than  by  Mc.  For  large  /?,  greater  than  30 
deg.,  this  could  bring  drastic  changes  to  the  flowfield 
since  a  structure  with  Me  >  1  may  effectively  be  sub¬ 
sonic.  Morkovin11  proposes  that  oblique  structures  be¬ 
gin  dominating  the  flow  field  once  Mc  becomes  super¬ 
sonic. 

The  best  way  to  to  visually  determine  flow  three- 
dimensionality  is  to  obtain  a  plan  view  of  the  shear 
layer,  which  necessitates  windows  on  the  top  and  bot- 
ton  walls  of  the  test  section.  In  the  current  facility, 
such  windows  would  have  been  difficult  to  install.  An 
easier  approach,  though  much  inferior,  is  to  visual¬ 
ize  the  flow  in  a  tilted  direction  through  the  existing 
sidewall  windows.  This  approach  was  tried  and  is  de¬ 
picted  on  Fig.  10.  The  apparatus  was  rotated  about 
its  flow  axis  so  that  the  schlieren  beam  entered  the  test 
section  at  a  12-deg.  inclination.  Due  to  the  resulting 
optical  distortion,  structure  obliquity  was  greatly  am¬ 
plified  on  the  image.  For  example,  a  30-deg.  obliquity 
would  have  appeared  as  70-deg.  obliquity  on  the  pho¬ 
tos.  Therefore,  only  obliquity  that  exceeded  70  deg. 
on  the  photo  could  be  considered  substantial  in  accor¬ 
dance  with  (5). 


Tilted-view  experiments  were  tried  with  cases 
H17N28  and  H26N28.  The  knife  edge  was  oriented 
perpendicular  to  the  streamwise  direction  to  empha¬ 
size  gradients  in  that  same  direction.  Although  the 
pictures  suffer  from  lack  of  clarity,  inherent  to  the  high 
distortion,  features  of  structures  were  picked  up.  By 
following  their  evolution  through  the  two-spark  sys¬ 
tem,  such  features  were  distinguished  from  stationary 
shock  waves  that  would  otherwise  have  created  confu¬ 
sion. 

The  structures  do  not  appear  to  be  convecting  in 
a  highly-oblique  fashion.  Obliquity  on  the  photos  did 
not  exceed  45  deg.,  corresponding  to  15-deg.  real  obliq¬ 
uity.  Fig.  11  shows  a  two-spark  photo  of  case  H17N28. 
The  trailing  edge  is  at  the  left  boundary.  Arrows  indi¬ 
cate  a  spanwise  structure  and  its  later  position.  That 
structure  looks  roughly  two-dimensional  and  is  repre¬ 
sentative  of  structures  seen  in  other  photos  of  that  case 
and  in  photos  of  case  H26N28.  The  convective  veloc¬ 
ity  computed  from  the  photo  agrees  with  that  obtained 
from  the  untilted  experiments  (Fig.  5). 

It  must  be  emphasized  that  visualizing  the  flow  with 
such  large  optical  distortion  is  not  a  reliable  means  of 
detecting  three  dimensionality.  The  results  presented 
in  this  section  should  therefore  be  treated  as  very  pre¬ 
liminary. 


Trailing-Edge  Devices 

Equation  (5)  indicates  that  if  the  flow  is  made  highly 
three  dimensional,  its  convective  Mach  numbers  would 
effectively  be  reduced.  From  Fig.  2,  it  is  evident  that 
this  could  cause  the  growth  rate  to  increase.  Sev¬ 
eral  experiments  were  conducted  with  the  the  goal 
of  enhancing  three  dimensionality,  and  hence  mixing, 
through  trailing-edge  modifications.  Three  types  of 
trailing-edge  devices  were  tried:  a  vortex  generator; 
trip  wires  oriented  at  60  deg.  to  the  trailing  edge;  a 
sawtooth  extension  of  the  trailing  edge,  with  teeth  in¬ 
clined  at  40  deg.  They  are  shown  schematically  in  Fig. 
12  with  relevant  dimensions. 

The  devices  were  applied  to  case  H17N28.  The 
boundary-layer  momentum  thicknesses  immediately 
upstream  of  the  devices  were  calculated  to  be 
01=0. 1mm  and  02=0. 05mm,  approximately.  On  the 
growth-rate  versus  Mc  curve  (Fig.  2),  this  case  is  lo¬ 
cated  at  the  kink  where  the  growth  rate  bottoms  out. 
Thus  a  decrease  in  Me<//,  even  a  small  one,  should 
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have  increased  the  growth  rate  by  a  noticeable  amount. 

None  of  the  devices  enhanced  the  growth  rate  by 
more  than  5%.  To  illustrate  the  insensitivity  of  the 
flow  to  the  mixing-enhancement  attempts,  Fig.  13 
shows  a  photo  of  case  H17N28  with  and  without  vortex 
generator.  The  growth  rates  appear  almost  identical. 
Tilted  views  of  that  flow  did  not  appear  different  from 
the  cases  with  clean  trailing  edge. 

Strange  Waves 

As  mentioned  earlier,  schlieren  photos  with  the  knife 
edge  perpendicular  to  the  flow  direction  were  also  pro¬ 
duced  in  the  experiment.  This  knife-edge  orientation 
picks  up  gradients  in  the  streamwise  direction,  which 
are  generally  much  weaker  than  those  in  the  transverse 
direction.  It  gives  a  more  detailed  view  of  turbulence, 
since  the  bias  from  the  refractive-index  difference  be¬ 
tween  the  two  streams  is  effectively  removed.  In  fact, 
it  shows  that  the  shear  layer  becomes  turbulent  im¬ 
mediately  downstream  of  the  trailing  edge  (Fig.  14), 
something  that  is  hard  to  establish  from  pictures  taken 
with  a  parallel  knife  edge  (Fig.  13).  It  was  mainly  ap¬ 
plied  to  cases  H17N28  and  H26N28. 

A  peculiar  phenomenon  was  observed  when  visual¬ 
izing  case  H26N28  in  that  fashion.  Waves  that  were 
roughly  normal  to  the  flow  direction  and  of  smooth 
texture  appeared  consistently  on  every  photo  of  that 
case.  Fig.  14  shows  a  typical  one.  The  waves  extend 
into  both  free  streams  and  are  relatively  stationary. 
They  did  not  appear  in  case  H17N28,  perhaps  because 
they  did  not  exist  or  because  they  were  too  weak. 

Although  these  waves  cannot  yet  be  fully  explained, 
it  is  useful  to  examine  what  they  are  not:  they  are  not 
normal  shock  waves,  or  the  pressure  rise  across  them 
would  have  been  enormous;  they  are  not  radiating  di¬ 
rectly  from  the  structure,  or  a  regular  pattern  with 
slope  equal  to  the  Mach  angle  at  MCl  =  1.47  (43  deg.  ) 
would  have  formed.  They  are  not  caused  by  gradi¬ 
ents  external  to  the  test  section,  because  photos  with 
no  flow  during  the  same  experiment  showed  a  uniform 
background. 

These  waves  are  thought  here  to  be  caused  by  an  in¬ 
tricate  interaction  between  a  supersonic  structure  and 
the  test-section  walls  and  their  boundary  layers.  The 
fact  that  they  appear  in  case  H26N28  with  A/c,  =  1.47 
and  they  do  not  appear,  or  are  too  weak,  in  case 
H17N28  with  A/Cl=0.83  supports  that  hypothesis.  A 


structure  with  Me,  =  1.47  emits  Mach  waves  that  reach 
the  walls  and  reflect  back  on  the  main  flow,  creat¬ 
ing  a  complex  pressure-wave  system.  A  structure  with 
AfCl=0.83  cannot  emit  Mach  waves  except  within  its 
close  proximity  where  the  Mach  number  is  supersonic. 

If  the  waves  are  indeed  generated  in  this  manner, 
it  is  still  unclear  why  they  do  not  convect  with  the 
structure  but  remain  almost  stationary.  It  is  possible 
that  they  are  initially  created  by  a  mechanism  similar 
to  the  one  described  above,  but  then  acquire  a  nature 
of  their  own.  Obviously,  more  experiments  are  needed 
to  explain  this  strange  phenomenon. 

Discussion 

The  experimental  results  raise  the  question  as  to  what 
the  appropriate  compressibility-effect  parameter  for 
the  shear  layer  is.  Previously,  the  isentropic  model 
with  MCl  «  MC3  made  the  answer  easy:  either  Mc  was 
a  suitable  choice.  Now  that  that  model  is  found  to  be 
inaccurate,  that  choice  must  be  reexamined.  The  re¬ 
sults  of  this  study  suggest  two  ways  to  describe  shear- 
flow  compressibility  in  a  consistent  manner.  One  way 
is  to  use  the  larger  of  the  two  Mc's.  Since 

the  lower  Mc  is  found  here  to  be  always  low  subsonic, 
its  contribution  is  secondary.  Thus,  MCm„  directly  de¬ 
fines  the  compressibility  associated  with  the  large-scale 
structure.  An  estimate  for  Mc„„  can  be  obtained 
from  Fig.  8,  assuming  that  the  curve  has  universal 
character.  That  leads  to  the  other,  more  indirect,  way 
of  describing  compressibility  by  the  parameter  -. 
As  already  noted,  it  determines  when  a  structure  has 
to  be  supersonic  regardless  of  the  value  of  Uc.  Its  good 
correlation  with  MCm„  suggests  that  it  too  is  a  consis¬ 
tent  measure  of  compressibility.  A  plot  of  growth  rate 
versus  -A^-  would  be  virtually  indistinguishable  from 
the  plot  of  Fig.  2.  Plotted  versus  MCm„,  the  growth 
rate  would  still  fall  on  a  universal  curve,  provided  that 
the  relation  of  Fig.  8  is  itself  universal. 

The  two-spark  schlieren  photos  revealed  some  re¬ 
markable  features  of  the  compressible  shear  layer. 
Generally,  the  large-scale  structure  appears  frozen 
from  one  exposure  to  the  next,  even  though  the  ve¬ 
locity  difference  across  the  structure  is  of  the  order 
of  the  convective  velocity  and  the  structure  has  typ¬ 
ically  travelled  two  of  its  body  lengths.  There  is  no 
evidence  of  pairing  or  coalescence.  It  is  thus  of  great 
interest  to  inquire  where  these  structures  are  “born”. 
If  they  are  generated  near  the  trailing  edge,  and  then 
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merely  convected  and  expanded  downstream,  the  flow 
field  in  the  vicinity  of  the  trailing  edge,  which  in¬ 
cludes  the  wake  region,  may  influence  their  behav¬ 
ior.  The  difference  between  supersonic-supersonic  and 
supersonic-subsonic  combinations,  noted  above,  could 
be  connected  to  the  fact  that  the  wake  zone  is  more 
dominant  in  the  former  cases  than  in  the  latter  ones. 

Although  the  revised  structure  model  makes  use 
of  shock  waves  emanating  from  the  structures,  such 
waves  are  not  clearly  seen  in  the  schlieren  photos.  This 
may  be  due  to  the  rapid  weakening  of  their  strength 
and  their  transformation  to  weak  waves  by  action  of 
the  surrounding  expansion  waves  (Fig.  7).  In  fact,  if 
the  shock  waves  retained  their  strength  far  from  the 
structure,  they  would  greatly  alter  the  free  stream  con¬ 
ditions:  U\  would  have  to  be  continuously  decreasing 
and  U2  continuously  increasing,  leading  to  an  impos¬ 
sible  situation  for  a  shear  layer.  In  the  experiment, 
(/ 1 ,  U2,  and  Uc  were  fairly  independent  of  streamwise 
location.  From  that  point  of  view,  it  is  not  surprising 
that  such  shock  waves  do  not  extend  far  into  the  free 
stream. 

The  asymmetric  trends  observed  in  the  experiments 
could  have  a  significant  impact  on  mixing  and  com¬ 
bustion.  The  mass  entrainment  ratio 

£  _  £2  Uc  -  Ui  _  IP2  72 

Pi  U\  Ue  ~~  y  pi  7!  MCi 

is  far  from  unity  in  the  highly-compressible  cases  due 
to  the  large  difference  between  MCl  and  AffJ.  The 
uncertainty  in  predicting  which  way  the  asymmetry 
will  occur,  hence  which  side  will  be  preferentially  en¬ 
trained,  may  pose  challenges  for  the  design  of  efficient 
supersonic  combustors. 

These  asymmetries  are  not  without  theoretical 
precedent.  Similar  trends  are  found  in  the  inviscid 
linearized  stability  analysis  for  an  unbounded  shear 
layer  with  uniform  density  by  Blumen  et  a/.12:  as  the 
convective  Mach  number  exceeds  1,  the  disturbance 
phase  velocity  departs  from  the  mean  velocity  and  ap¬ 
proaches  either  U\  or  C/2  -  Presence  of  solid  boundaries 
in  the  flow  can  introduce  new  modes  of  instability, 
some  of  which  are  asymmetric.  Such  modes  are  ev¬ 
ident  in  Tam  and  Hu’s  13  eigenfunction  distributions 
for  a  bounded  shear  layer. 

Given  the  radiative  nature  of  a  supersonic  distur¬ 
bance,  the  effects  of  the  surrounding  walls  cannot  be 
neglected,  as  is  frequently  done  in  the  subsonic  case. 
The  strange  waves  seen  in  the  schlieren  photos  of  Fig. 


14  are  perhaps  a  hint  that  walls  have  to  be  an  integral 
part  of  any  realistic  compressible  shear  layer  model. 

The  failure  of  mixing  enhancement  through  efforts 
to  increase  flow  three-dimensionality  may  have  two  dif¬ 
ferent  and  highly-speculative  interpretations:  (a)  the 
flow  instability  is  basically  two-dimensional  and  resists 
any  departure  from  that  state;  (b)  the  instability  is 
highly  three-dimensional  to  begin  with,  in  which  case 
the  devices  would  have  little  to  contribute.  The  vi¬ 
sualization  of  distinct  large-scale  structures  through 
the  schlieren  system,  which  gives  a  span  wise-integrated 
view  of  the  flow,  and,  to  a  lesser  extent,  the  tilted-view 
experiments,  tend  to  support  argument  (a).  However, 
this  issue  cannot  be  settled  until  a  direct  plan  view  of 
the  compressible  turbulent  shear  layer  is  obtained.  It 
must  also  be- noted  that  this  study  has  by  no  means 
exhausted  all  the  methods  for  efficient  mixing  enhance¬ 
ment.  Other  configurations  may  hold  more  promise. 


Conclusion 


Inclusion  of  shock  waves  in  the  model  for  a  com¬ 
pressible  large-scale  structure  gives  rise  to  highly- 
asymmetric  solutions  for  the  convective  Mach  num¬ 
bers.  This  helps  interpret  the  large  difference  between 
Mc,  and  Mc,  seen  in  the  experiments.  The  parameter 
af+Zi  delineates  the  likelihood  of  asymmetric  solutions 
ana  is  a  broad  measure  of  shear-flow  compressibility. 
Tilted-view  experiments,  although  of  preliminary  na¬ 
ture,  did  not  reveal  pronounced  obliquity  of  the  large- 
scale  structure.  The  failure  of  mixing-enhancement 
attempts  through  trailing-edge  devices  underlines  the 
inherent  stability  of  highly-compressible  shear  layers. 
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Fig.  1  Shear-layer  geometry  and  nonmenclature. 
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Fig.  2  Isentropic  model  of  large-scale  structure. 


Fig.  3  Normalized  growth  rate  versus  theoretical 
M'X  (Refs  2,9). 


Fig.  4  Operating  principle  of  two-spark  schliercn 
system. 
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Fig.  5  Two-spark  schlieren  photos,  covering 
flow  from  50mm  to  140mm  downstream  of  trailing 
edge, 
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Fig.  8  Dependence  of  MCl  and  Mc,  on  Uc  for 
arbitrary  conditions.  M *  indicates  the  level  to  be 
exceeded  for  asymmetric  situations  to  be  likely. 


Fig.  6  Experimental  and  theoretical  Me,  —  vs.  — 
MCj  plot  for  the  ten  cases  investigated  here. 
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Fig.  7  Revised  compressible  large-scale-structure 
model.  Sketch  of  total-pressure  distribution  ( Pt ) 
along  streamlines  leading  to  stagnation  point  (0) 
illustrates  the  symmetric  and  asymmetric  possi¬ 
bilities. 


Fig.  9  The  maximum  of  experimental  MCl  and 
Mc2,  MCm„,  versus 


Fig.  10  Tilted- view  method  for  visualizing  span- 
wise  structure.  Flow  is  coming  out  of  the  paper. 
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1I17N28:  (a)  with  vortex  generator;  (b)  without. 
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Fig.  11  Tilted  view  of  case  H 17N28  with  Uc  mea¬ 
surement. 


Fig.  12  Trailing-edge  modifications:  (a)  vortex 
generator;  (b)  slanted  trip  wires;  (c)  saw-tooth 
extension. 
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Fig  M  Two-spark  photos  of  case  H2fiN28  taken 
with  knife  edge  vertical.  Upper  photo  covers 
near  field  (0  to  SHInim),  where  f,  measurement 
is  shown;  lower  photo  covers  far  field  (50mm  to 
MOmm).  Strange  waves  are  visible  in  both  pho¬ 
tos. 


